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Abstract 

The recently found non-BPS multi-wall configurations in the M — \ supergravity in four 
dimensions is shown to have no tachyonic scalar fluctuations without additional stabilization 
mechanisms. Mass of radion (lightest massive fluctuation) is found to be proportional to 
^g-TrAi?/2^ where A is the inverse width of the wall and R is the radius of compactified 
dimension. We obtain localized massless graviton and gravitino forming a superniultiplet 
with respect to the Killing spinor. The relation between the bulk energy density and the 
boundary energy density (cosmological constants) is an automatic consequence of the field 
equation and Einstein equation. In the limit of vanishing gravitational coupling, the Nambu- 
Goldstone modes are reproduced. 
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1 Introduction 



In the brane- world scenario 0E1E1, our four-dimensional world is to be realized on topological defects 
such as walls. To obtain realistic unified theories beyond the standard model, supersymmetry (SUSY) 
has been most useful Moreover, SUSY helps to construct topological defects like walls as BPS 
states [H] that preserve part of SUSY. For a reahstic model, understanding SUSY breaking has been an 
important problem, which is addressed in the SUSY brane- world scenario extensively 0-^31 ■ Models 
have been constructed that realize one such idea : coexistence of BPS and anti-BPS walls produces SUSY 
breaking automatically ^J). In particular, the SUSY breaking effects are suppressed exponentially as 
a function of distance between walls. On the other hand, non-BPS multi-wall configurations are not 
protected by SUSY and need not be stable. Such non-BPS wall configurations was successfully stabilized 
by introducing topological quantum numbers, such as a winding number |15[ I16| . The physical reason 
behind the stability is simple : a BPS wall and an anti-BPS wall with winding numbers generally exert 
repulsion, which then pushes each other at anti-podal points of the compactified dimension. 

One of the most attractive models in the brane-world scenario is the model with the warped metric 
[51 |S| . A possible solution of the gauge hierarchy problem was proposed in the two brane-model [2 , and 
a localization of graviton on a single brane was found even in a noncompact space f3| at the cost of fine- 
tuning between bulk cosmological constant and boundary cosmological constant at orbifold fixed points. 
Supersymmetrization of the thin- wall model has also been constructed in five dimensions jl7|-|19j. It 
is natural to ask if the infinitely thin branes in these models can be replaced by physical smooth wall 
configurations made out of scalar fields |2ni~Ell- We have succeeded in constructing BPS as well as 
non-BPS solutions in the M = 1 supergravity coupled with a chiral scalar multiplet in four dimensions 
123] ■ A similar BPS solution has also been constructed in five-dimensional supergravity |25l I2f)j . In the 
limit of vanishing gravitational coupling k — > 0, our model reduces to the model having the exact solution 
of non-BPS multi- walls JSj. Therefore the model is likely to be stable thanks to the winding number 
near the weak gravity limit. However, we need to address the issue of stability in the presence of gravity, 
since the radius of the extra dimension is now a dynamical variable which might introduce instability 
into the model. There have been a number of works to analyze the stability of the infinitely thin wall 
[2H1-[12, especially in the presence of a stabilizing mechanism due to Goldberger and Wise [33) . 

The purpose of our paper is to study the stability of the model with winding number in the presence 
of gravity and to analyze the mass spectrum of fluctuations on the BPS and non-BPS solutions. We find 
that there are zero modes of transverse traceless fluctuations localized on the wall which play the role of 
the graviton in our world on the wall. The BPS solution has also gravitino zero mode which is localized 
on the wall and forms a supermultiplet with the graviton under the surviving supergravity transformation 
with the Killing spinor of the BPS solution. We obtain that the BPS solution has no other zero modes, 
and no tachyonic fluctuations. For instance, we find that possible additional massless tensor and scalar 
modes are either gauge degrees of freedom or unphysical (the mode function is not normalizable) . As for 
the non-BPS solution, we find that another possible zero modes of the transverse traceless fluctuations 
of metric can be gauged away and that there exists no zero mode other than the graviton localized on 
the wall. To obtain a concrete estimate of the mass spectrum, we need to use approximations. We 
use small width approximation where the width A^^ of the wall is small compared to the radius R 
of compactified extra dimension. We find that the non-BPS solution has no tachyonic fluctuations in 
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spite of the dynamical role played by the radius of the compactified dimension. Tensor as well as scalar 
fluctuations have massive modes, without any tachyons. This result shows that our non-BPS solution is 
stable without introducing an additional stabilizing mechanism such as the Goldberger-Wise mechanism 

The lightest massive scalar mode is usually called radion. We can evaluate the mass of the radion on 
our non-BPS background at least for R 3> A^^, where R is the radius of the compactified dimension and 
is the width of the wall. We find that the mass squared of the radion is given by 

ml cx A^e-^^^ (1.1) 

It is interesting to note that the mass scale is given by the inverse wall width A, and that it becomes 
exponentially light as a function of the distance ttR between the two walls. This behavior is precisely the 
same as the previous model in the global SUSY case [TS] . 

Modes of fermions including gravitino are also analyzed. We find that the Nambu-Goldstone modes 
can be reproduced in the limit of vanishing gravitational coupling both for bosonic and fermionic modes. 

Our BPS solution has a smooth limit of thin walls where it reproduces the Randall-Sundrum model 
[24|. In the original Randall-Sundrum model, the fine-tuning was necessary between the boundary and the 
bulk cosmological constants. However, the necessary relation between bulk and boundary cosmological 
constants is now an automatic consequence of the equation of motion of scalar fields and Einstein equation 
in our model. We no longer need to impose a fine-tuning on input parameters of the model. 

Sec. 2 summarizes our model and solutions briefly. Sec. 3 separates various bosonic modes with respect 
to the surviving Lorentz symmetry (tensor and scalar modes) and addresses the question of stability of 
the BPS solution. Sec. 4 discusses the stability of non-BPS solution and evaluates the mass of the radion. 
Sec. 5 deals with the fermionic modes. The gauge fixing to the Newton gauge is justified in Appendix A, 
and some illustrative cases of potential in the conformal coordinate are worked out in Appendix B. 

2 Brief review of BPS domain wall in SUGRA 

2.1 Lagrangian and BPS equations 

We consider a chiral multiplet containing scalar and fermion x with the minimal kinetic term and 
the superpotential P, and the gravity multiplet containing vielbein Cm- and gravitino V'm"- The local 
Lorentz vector indices are denoted by letters with the underline as a, and the vector indices transforming 
under general coordinate transformations are denoted by Latin letters as to, 7i = 0, . . . , 3. The left(right)- 
handed spinor indices^ are denoted by undotted (dotted) indices as a{a). Then the Af = I supergravity 

^We follow conventions of Ref . 1341 for the spinor and other notations. 
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Lagrangian is given in four-dimensional spacetime as |84j 



2k2 

^2 _ _ 



-e 2 



(p*V„a™V„ + PV;™a'""^„) 



(2.1) 



where the gravitational coupling k is the inverse of the four-dimensional Planck mass Afpi, gmn is the 
metric of the spacetime and e is the determinant of the vierbein Cm— • The generalized supergravity 
covariant derivatives are defined as follows : 



T^rnX = drrtX + X^m ^Im[(/)*(9,„0] X, 

IK^ 

~2 



i'mi'n = dmll^n + V'n^m + ^Im[(/)*9„^] V'n, (2.2) 



where (jj„i is the spin connection and we use the notation Im[X] = in what follows. The scalar 

2i 

potential in the supergravity Lagrangian (|2.1() is given by 

V{4>, r) = e"'^*^ {\D^P\^ - 3«;2|F|2) . (2.3) 
The above Lagrangian (|2.1(l is invariant under the supergravity transformation : 

• 9 

8c_1p,n = 2K-^VrnC + i Kc"^ '''' P (J ,nC, " yCrm«CXCr"X " —'^M(t>*k^] '<Prn, 

■ 2 

Sex = iV2a"XVm^-V2e4r<^D^,P*C+^-lm[rSc4>]x, 
where ^ is a local SUSY transformation parameter and the covariant derivatives are given by 

^'m'/' = d,n(f> - -^K'lpmX, 
^ . 2 

2?mC = dmC, + C^m + — Im[(/)*a,„(/)] 



(2.4) 



(2.5) 



Next we turn to derive the equations of motion for solutions which depend on only one "extra" 
coordinate — y under the warped metric Ansatz 

ds^ = .9„„dx'"da;" = e^'^^y^'ij^^dx^dx" + dy"^ (n, iy = Q,l, 3), (2.6) 
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where Greek indices fi = 0, 1,3 denote three-dimensional vector transforming under general coordinate 
transformations, and rj^i, — diag(— ,+,+) denotes three dimensional flat spacetime metric. All the 
geometrical quantities can be written in terms of the function A{y) in the warp factor and its derivatives 
with respect to the extra coordinate y. For later convenience, we write formulas in general D space-time 
dimensions in the following : 

1. vierbein 

e™^ = diag (e^, e^^, 1, e^) , ej^ = diag (e^^, e"^, 1, e"^) , (2.7) 

2. spin connection 

{X^m)a = ^t^mob ^mab = A (Sa^ebm ~ Sb^eam) , (2.8) 

3. Ricci tensor 

Rrnn = e^-^ (^A + {D - 1)A^^ Vt^.^m'^Sn'^ + (D - 1) (^A + i^) 6r„^Sj, (2.9) 

where a dot denotes a derivative with respect to y, A = dA/dy, and we turn off all the fermionic fields 
as a tree level solution. The energy momentum tensor is given in terms of the scalar potential V{(j3, (jf) 
in (ESI 

Tmn = d,n4>*dn4> + d^<t>dn(t>* ~ 5,„n {g""' dk4>* di<t> + V{<t>, 4>*)) ■ (2.10) 
Plugging these into the Einstein equation^ Rmn = ^i^'^Tmn, we obtain 

^^- (^ 2) (2-11) 

The field equation for the scalar (j) in the chiral multiplet takes the form : 

■ci,+ {D~l)A4>=^. (2.12) 

Notice that only two out of the three equations in Eas. H2.11() and l|2.12|l are independent (assuming only 
one real component, say the real part of the scalar field <j) is nontrivial in the solution). Any one of three 
equations are automatically satisfied if others are satisfied. 

It is well known that special type of solutions for these nonlinear second order differential equations 
are obtained as solutions of a set of the first order differential equations, the so-called BPS equations 
which guarantees the partial conservation of SUSY. Similarly to the global SUSY case, the BPS equations 
can be derived from the half SUSY condition where we parametrize the conserved SUSY parameter as 

ay)=e'^y^aH{y). (2.13) 

That is, we demand that the bosonic configuration should satisfy 5qx — ^Q^'m — for the parameter C(y) 
in Ea. (|2.13|l . The BPS equations for the metric are derived from the condition for the gravitino. From 
TO = /I = 0, 1, 3 components the first order equation for the warp factor A is derived : 

- = «"^e^ [e**^ + iK^e'^'''''l'p\ ct^C, (2.14) 



2We define fWi„ = Tm„ - -fT=^g^-aT^ k- 
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A 



(2.15) 



From m = 2 component we find the first order equation for the Killing spinor Q 



= 5qiIj2 = 2k' 



c 



-Im 



c- 



IK 



"e 2 



(2.16) 



Rewriting the half SUSY condition as Cq = e^C^+s)]^^] and substituting it into the above equation, we 
find 123 



i4i = |iy, 



-K^Im 



(2.17) 



On the other hand, the first order equation for the matter field (p is derived from the half SUSY condition 
for the matter fermion x • 



-ze e 2 



(2.18) 



Ea. (|2.15|l . H2.17|l and H2.18(l are collectively called BPS equations. One can easily show that solutions of 
the BPS equations satisfy the equations of motion H2.11|l and (I2.12|l . Notice that the Ea. H2.18|l and the 
second equation of Ea. (|2.17|l do not contain the metric, so we can solve this as if the scalar field decouples 
from gravity. Once the configuration of the scalar field 4> and the phase 9 are determined, the warp factor 
A is obtained from Ea. H2.15|l . Finally, the Killing spinor ^ is also determined from the first equation of 
Eq.(jnZ|)- 



2.2 Exact BPS solution 

Recently, we found the exact BPS solutions for the periodic model in SUGRA |24j, by allowing the 
gravitational correction for the superpotential as follows 



P(0) = e""^ 



A3 



9 



9 



- sm ■ ^ 

2 A 



(2.19) 



where A is a coupling with unit mass dimension and g is a dimensionless coupling. We introduced this 
modification for the superpotential in SUGRA to maintain the periodicity of the model with the aid 
of the Kahler transformation. This modification for the superpotential gives SUSY vacua which do not 
depend on the gravitational coupling n. This was crucial for us to obtain the exact BPS solutions in 
SUGRA. The superpotential (|2.19() yields the following scalar potential : 



V 



A4 



9 

a' 



Im (7> sm— f 

9 ' ' A 



3k2A2 



9 

a' 



(2.20) 



The SUSY vacua are determined from the condition D^P = 0. For the above modified superpotential 
we find that the SUSY vacua are periodically distributed at = — ^— + m:^ , (n € Z) on the real axis 
in the complex (jj plane. 

In order to determine the scalar field configuration, we need to solve the second equation of Eq. (|2.17(l 
for 6 together with the equation for scalar field : 

'— cos -7-1; 
9 [ A 



2iAK2 g 

Im <P sm 

9 A 



(2.21) 
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To solve Eas. (|2.17|l and (|2.21|l . we choose = Ini[(/)] — and = ±— at a point, say y = yi as an initial 
condition for the imaginary part (f>i{y) of the scalar field and the phase 0{y). Then these equations tell 
that </>! = 9 = at y — yi. Therefore we find 0i = and 9 — zL— at any y. At this stage, only the 
real part of (p has a nontrivial configuration in the extra dimension y. We shall call those scalar fields 
that have nontrivial configurations as a function of the coordinate of extra dimension, as "active" scalar 
fields. The scalar potential along (pi ^ surface is given by the following potential Vr for the real part 
(pYi = ~^ of the scalar field : 



A" 

^r(0r) = — 
9^ 



7 9 , 
cos-" ^<?!)R, - 



3k2A2 



•25, 

sm — 



(2.22) 



It has been shown that the following form of scalar potential with a real "superpotential" P((/)r) of a real 
scalar field 0r ensures the existence of a stable AdS vacuum in gravity theories in D dimensions |35[ 



D 



D 



(2.23) 



if there is a critical point in P^cpn), even though supersymmetry is not required in this form. Let us 
note that our scalar potential is compatible with the above form of the scalar potential. In our case, the 
"superpotential" is given by 



A3 g 
9 A 



A /TT 



mr] , our scalar potential 1^,(0r) in Eq. 



Since this P has critical points at (pfi 
critical points as stable AdS vacua. 

The remaining BPS equations for the active scalar field and the warp factor are of the form 



(2.24) 
has these 



i , dP A' g ^ 

<)R = ± V— = ± cos -0R, 

d0R g A 



A = tk'^P = T 



k2A3 



9 



9 



2 sm-(?;R. 



(2.25) 
as an initial 



Let us solve these BPS equations by choosing a SUSY vacuum 0r = — ('ip(— 1)"^ + mr 

9 ^ 2 

condition at y = — oo. We shall consider the solution for the BPS equations H2.25|l with the sign correlated 
to the sign of the initial condition at y = — oo. The exact BPS solutions are found to be of the form : 
A 



OR 



9 



(-!)"{ 



2tan-ie=^^(^-^''' 



e = [coshA(?; - yo)]" 



(2.26) 



where k 



9^ 



lutions interpolate between the two SUSY vacua, from 0r = — ^t(~1)"'^ + ?^7r^ at y = — oo to 



— — is the inverse of the curvature radius of the AdS spacetime at infinity. These so- 

A 

^ 9 
0R — ^±(—1)"— + rnr^ at y — +oo. We denote yo the modulus parameter of these solutions and we 
suppress an integration constant for A which amounts to an irrelevant normalization constant of metric. 
Eq. (|2.17|) determines the Killing spinors which has two real Grassmann parameters ei , £2 corresponding 
to the two conserved SUSY directions on the BPS solution'^ : 



c 



C 2 X 



£2 



(2.27) 



These Killing spinors are the corrected results of those in our previous work [an. 
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Ub+^\ A i i[coshA{y-yo)] , for = |, 
e2(«+2 jg- = ,^ _^ 2^ (2.28) 

[coshA(2/-yo)] , for 61 = --. 

Our model has a smooth hmit of thin waUs where it reproduces the Randall-Sundrum model |24j . 
Notice that we do not need any fine-tuning of input parameters of the model, in contrast to the original 
Randall-Sundrum model. The necessary fine-tuning between bulk and boundary cosmological constants 
is now an automatic consequence of the equation of motion of scalar fields and Einstein equation in our 
model. 

2.3 non-BPS solution 

Assuming that only single real scalar field has nontrivial classical configuration, the equations (|2.11(l 
and (|2.12|l reduce to 

2^2 .„ .„ 2^2 / :o \ ■; . : ^ f dVr 



{^l-V^), 4>^ + {D-l}A,Pj, = -^. (2.29) 



D-2^^' {D-l)iD-2) V'^ "7' ^ ' 2d(t) 

It has been shown that the above set of coupled second order differential equations is equivalent to the 
following set of nonlinear differential equations |22l I28| . Given the scalar potential Vfi(0r,), we should 
find a real function VF(0r) by solving the following first order nonlinear differential equation 



dWi^bn) ,2 L, {D-l){D-2) 

^ ^^^V^'^^^^ + 2 -'W'iM- (2.30) 

Then 0R(y) and A{y) are obtained by solving the following two first order differential equations 

^-(^) = V^^' My) = -^'wiM- (2.31) 

If we choose the "superpotential" P as a real function W, H2.30|l and (|2.31|) are satisfied by the scalar 
potential H2.23|l and the BPS equations H2.25|l . Therefore these set of first order nonlinear differential 
equations includes all the BPS solutions as part of the solutions. However, it is important to realize that 
()2.30(l and H2.31|l are equivalent to the set of Einstein equation and the scalar filed equation, and hence 
give all the non-BPS solutions as well. 

We have been able to construct non-BPS multi-wall solutions to the Einstein equation (|2.11|l and the 
field equation (|2.12() using the above method of nonlinear equations i24|. We have also found that BPS 
solutions are the only solution that do not encounter singularities at any finite y. To obtain any other 
regular solution, especially non-BPS solutions, negative cosmological constant has to be introduced at 
some boundary. Since we are interested in periodic array of walls where extra dimension can be identified 
as a torus with possible division by discrete groups (orbifolds), we introduced the cosmological constant 
and obtained a number of interesting non-BPS solutions |24j . 

The above nonlinear differential equation (|2.30(l gives a set of solution curves which fill once and only 
once the entire (0r, W) plane except forbidden regions defined by Vr + {D — 1){D — 2)k'W'' /2 < 0. Let 
us denote the solution curve starting from an initial condition Wq at (/)r — (/jr^q as VF((/)r; (0r,o, W^o)). 
A boundary cosmological constant Xi at yt gives a jump of derivative of the function A{y) in the warp 
factor. Let us denote the value of the scalar field at the boundary yi as 0R,i. Ea. H2.31|) shows that this 
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jump of A{y) is satisfied by cutting the solution curve and jump to another solution curve at (/)p{_,; with 
the constraint 

A, = 2 {Wi<j)n,^ + Wi^j^n., - e)) • (2.32) 

Since we are interested in minimum amount of inputs at boundaries, we wish to implement only the 
boundary cosmological constant without any boundary potential for scalar fields 0fi, contrary to many 
other approaches characteristic of the Goldberger-Wise type of the stabilization mechanism [331, |22| . 
|23|. Therefore we need to maintain the derivative dW/dcf)^ to be smoothly connected at the boundary. 

Since Ea. H2.3()|l gives the same value of derivative dW/dip-[{ for ztW, we can connect the solution curve 
at any value of (/)p{^i if we switch from a solution curve going through 4>-b.j to another one going through 
— W, Ea. H2.32() gives the necessary cosmological constant at this boundary as A = 4:W{(j)^^i). There 
may be other possibilities to connect the solution curves, but this is the simplest possibility that covers 
many interesting situations. 

To be definite, we shall consider walls that have simple symmetry property under the parity Z2 '■ 
(/)p{ — 0R. Let us start a solution curve going through 0r = 0, Wo > 0. Then the solution curve goes 
above the forbidden region. To obtain a non-BPS solution which is odd under the Z2 transformation, 
we place a boundary at 0r — with a positive cosmological constant by an amount Ao = AWq > 0. On 
the other hand, we can place a boundary at any > with a negative cosmological constant Ai = 
—4W{(j)B.,i, ('/'R = 0, Wo)). However, we can obtain a multi-wall solution that have simple transformation 
property under the Z2 by placing another boundary at integer multiple of (/)r — Kt: /{2g). 




g0/A 




Figure 1: Non-BPS solution with unit winding number with the /{Z2 x Z2) symmetry. Discontinuities in 
W{'j)) at g4>/^ = 0, TT (7r/2,37r/2) correspond to positive (negative) boundary cosmological constants, (a) (f)W 
plane, and (b)The function A[y) in the warp factor. 



If we place the first boundary at the vacuum point (/)r = A7r/(2f;), we obtain a simplest model in 
the sense that the energy density at the second boundary at 0r — K'K/{2g) is purely made of negative 
cosmological constant 

A = -4W(0R = A7r/(2g)). (2.33) 

The magnitude of this negative cosmological constant becomes the same as the total energy of the 
wall centered at 0r = in the limit of large separation of two boundaries. Since the solution admits 
S^/{Z2 X Z2) symmetry, we call the coordinate at the second boundary y — tiR/2. The behavior of 
this non-BPS solutions in the W,4> plane is illustrated in Fig^^a). The corresponding function A{y) 
in the warp factor is illustrated in Fig^^b), where one should note that A{y) is linear near the second 
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boundary ai y = ■nR/2, showing that only the boundary cosmological constant exists apart from the bulk 
cosmological constant there. 




Figure 2: Non-BPS solution with unit winding number with the /Z2 symmetry. Discontinuities in W{(j)) at 
g4>/A. = (tt) correspond to positive (negative) boundary cosmological constants, (a) ())W plane, and (b)The 
function A{y) in the warp factor. 

As another solution, we can place the second boundary at 4>b. — Att/^, where the active scalar field 
develops another wall configuration. In this case, the negative cosmological constant — 4VF((/)r = Air/g) 
placed at the second boundary has magnitude which becomes twice the total energy of the wall centered 
at = in the limit of large separation of two boundaries. The behavior of this non-BPS solutions in 
the W, (f) plane is illustrated in Fig[2Ia). The corresponding function A{y) in the warp factor is illustrated 
in Fig|2Ib), where one should note that the function A{y) has additional kink behavior deviating from the 
linear exponent near the second boundary at y = ■nR, showing that there is an additional smooth positive 
energy density centered around the boundary besides the negative boundary cosmological constant in 
contrast to the previous /{Z2 x Z2) example in FigQ] 

3 Bosonic Fluctuation and the BPS Solution 

A Bogomolo'nyi bound has been derived for the energy density of the BPS domain walls in A/" = 1 SUGR A 
in four-dimensional spacetime [201 ■ They used the generalized Israel- Nester-Witten tensor, which was 
originally applied to a simple proof of the positive ADM mass conjecture in general relativity. However, 
the ADM mass may not be well-defined for domain walls, since they are extended to infinity. Therefore it 
is presumably still useful to check that there is really stability of the fluctuation on our wall configuration 
even in the case of BPS solutions. We shall present a general formalism to analyze the modes and their 
stability, and then apply it to the fluctuations around the BPS background conflgurations in this section. 
The equations and procedures obtained in this section can also be used to the non-BPS background 
solutions with appropriate additional inputs, which is dealt with in Sec0] 

3.1 Mode equations for the bosonic sector 

We start with the metric perturbation in the Newton gauge [3D], [32] : 

ds' = e^^ (r/^, + hj^ + dx^'dx" + (1 - 2{D - 3)B)dy\ (3.1) 
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where /iJJ is transverse traceless ?7'"'/iJJ = 0, d^h'J^J — 0. Some details for the procedure of this gauge 
fixing are given in Appendix A. This gauge is useful since the linearized equations become very simple. 
The linearized Einstein equations in D space-time dimensions (D = 4 in our specific model) read : 

(^e-^^Dn-i + dl + (D- l)Ady) h^^ = 0, (3.2) 
Dd-1 + dl + {iD - h)Ady + 2[D - 3) (p - l)i2 + i)) S = (3.3) 

(a^ + (i?-3)i)s = --^^«2^RVJR, (3.4) 

where the first line comes from the traceless part of (/i, v) component of the linearized Einstein equations, 
the second line from the trace part of (^, v) and the last from (^, 2) component. The (2, 2) component 
of the linearized equation is not shown, since it can be derived from Eas. H3.2| l- (|3.4(l . The linearized field 
equations give : 

j2t 



-2A 



+ dl + iD- l)Ady - = -2iD - 2),PndyB - [D - 3)^B, (3.5) 

^ 2 / d0R 

Ui - 0, (3.6) 



,-^^Dn-i + dl + {D-l)Ady-l^ 



2 dcf)^ background , 

where (^r(i) denotes the real (imaginary) part of the fluctuation of the scalar field (j) = 0R + '/'R + i'^i 
around the background fleld configuration 0r . Notice that the solutions of the linearized Einstein equation 
automatically satisfy the linearized field equations for the active scalar field (pn. Therefore, the Eas. H3.2|l - 
constitute the full set of independent linearized equations for the fields hJ^J , B, (^r and (pi. 



3.2 Tensor perturbation : localized massless graviton 

First we show that the linearized equation for the transverse traceless mode (graviton) given in Eq. (|3.2|l 
can be brought into a Schrodinger form. It can again be rewritten into a form of the supersymmetric 
quantum mechanics (SQM) which ensures the stability of the system. For that purpose we change the 
coordinate y into the conformally flat coordinate z defined as 

dz = e-^^y^dy, ds^ = e^^^^'' {ri^.^dx^'dx" + dz^) . (3.7) 

We also redefine the field as /iJJ = e^~"^/iJJ. In the following we use prime to denote a derivative in 
terms of z. Then the linearized equation H3.2(l becomes 

Uo-,hl^A^.^) = [-dl+V,{z)] hlJ{x,z), V,(z) = (^)'^'^ + (3.8) 

where Vt{z) is the potential in this "Schrodinger" type equation. For our BPS background solution H2.26(l 
the Schrodinger potential takes the form : 

fcA 



Vt{y) = [coshA(y-2/o)]" 



+ 2k'^ tanh^ A{y - yo) 



(3.9) 



cosh^ A{y-yo) 

where iT^ = 4g^^A^ is the tension of the wall and k = k?T^ . Although our model contains three 
parameters A, g and k, this potential depends on only two parameters k and A. If we take the thin wall 
limit where A — + cx3 fixing AT^ , we obtain (putting yo = 0) 

^ 25(z), tanh^ A?/(0) ^ 1, [cosh Ay(z)]~'^ ^ , A , 2 ' (3-10) 



cosh^Ay(z) {k\z\ + l) 
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with kz = sgn(y)e'^l^l — 1. Thus the Schrodinger potential (|3.9(l becomes precisely the potential of the 
Randall-Sundrum model : 

op 

Vt(^)^ , ~2k5{z). (3.11) 

{k\z\ + 1) 

We find that the part of action quadratic in h^"^ has no z dependent weight 

j dzd^'-'x v'^^v'^K^ (ao-i + dl - ivt^ h^;i, (3.12) 

in conformity with the absence of the linear term [221 i'^ i'^ ^^e Shrodinger type equation H3.8|l . We 
stress that this is written in terms of the conformal coordinate z and the redefined field . 

Defining mode equations by eigenvalue equations Htipniz) — [—9"^ + Vt(2:)] ipn{z) — m'^ipn{z) with 
mass squared eigenvalues m^, and assuming mode functions ipn{z) to form a complete set, the transverse 
traceless fields can be expanded into a set of effective fields /iJJ^"''(x) 

n 

Then the above quadratic action (|3.12(1 becomes 

S^Y.j ^n{z)Mz) ■ I d^'-'x r^V^^hU^^^ [On-i - ml] h'J^''\ (3.14) 

n.k 

Therefore the inner product for the mode function ^(z) should be defined as as 

(V'ilV'2) fdz Mz)Mz), (3.15) 



for which the usual intuition of quantum mechanics works. 

The Hamiltonian Ht can now be expressed in a SQM form as follows 

nt = QlQt, Qt = -d, + ^^-^A', Ql = d, + ^^^A', (3.16) 

where the "supercharge" Qt and QI are adjoint of each other at least for BPS background where no 
boundary condition has to be imposed. Therefore the Hamiltonian Tit is a nonnegative definite Hermi- 
tian operator"*, and its eigenvalues are nonnegative definite. Therefore we can conclude that the tensor 
perturbation has no tachyonic modes which destabilize the background field configurations at least for 
BPS solutions. 

There are two possible zero modes in the tensor perturbation. One is the state which is annihilated by 
Qtlhjj.J'"'^^ ) — 0, and another is the state defined as QI ^Qt|/iJJ' ') j = where Qt\h]j,J''-~^) ^ 0. Then 
zero modes are of the form : 

/.JJ(°)(^,z) = hl^(+Hx) e^-^W +/.;^J(-)(x) e^^(^) I dz e-(^-^)^(^), (3.17) 

where A(z) — A{y{z)). Notice that we must verify the normalizability of the wave-function to obtain a 
physical massless effective field in the case of noncompact space such as our BPS background. In the case 
of non-BPS background, the boundary condition has to be verified, which we shall consider in Sec^ The 

^Adjoint relation between Qt and qI and the Hermiticity of Tit are assured by the inner product defined in Eg. lij.l^ 
without z dependent weight. 
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first mode /iJJ^^'' in Ea. H3.17|l is normalizable if / dz e*^^^^^"^*^^^ < oo, corresponds to the graviton which 
is localized at the wall with a positive energy density, namely, if e^^~^^^ falls off faster than \z\~^ |29j. 
For our BPS solution H2.26|l the asymptotic behavior of the warp factor is of order \z\^^ . Therefore 
we obtain a normalizable massless transverse traceless mode h^i, which gives the physical graviton 
localized on the wall. 

On the other hand, the second term hjij^ ^ in Ea. (|3.17|l is not normalizable and is unphysical since 



ye^~ A J dz e^''^^^^^j ^ |z|'* at \z\ oo for our BPS solution. If there exists a regulator brane with 
a negative tension at some y, this mode can become normalizable and localizes at the negative tension 
brane in contrast to the graviton. If there are no bulk scalar fields (contrary to our model) as in the 
original Randall-Sundrum model of single wall, this zero mode corresponds to the physical massless field 
which was called radion in Ref.|27|. 

Our specific four-dimcnsional model of non-BPS wall gives a three-dimensional effective theory on 
the wall. Transverse traceless mode of graviton in three dimensions has no dynamical degree of freedom 
except possible topological modes. However, our formalism and analysis can be applied at each step to 
general I?-dimensional theories, once we obtain the relevant non-BPS solutions in such theories. In that 
respect, we believe that our findings should still be useful. 



-0 .1 



-0. 



1 




Figure 3: The Schrodinger potential for the tensor perturbation in terms of the coordinate system y. We choose 
yo = and k — 0.25. The solid line has A = 3 and the broken line has A = 1. 



The Schodinger potential can always be expressed in terms of y, but is difficult in terms of z explicitly^, 
since it is generally difficult to solve dz — e~^dy explicitly. If we express the potential in terms of y, we 
obtain a volcano type potential as shown in FiglHl The width of the well is ^ 2A^^ and the depth is 
~ kA. 

Next we turn to analysis of the massive Kaluza-Klein (KK) mode. There are no modes with negative 
mass squared in the tensor perturbation, as we have already shown. Since the Schrodinger potential 
()3.9|l vanishes asymptotically (z — ±oo), all the massive KK modes are continuum scattering states with 
eigenvalues 777? > 0. In order to examine the mode functions of the massive KK modes, we look into the 
region far from the wall, namely A|y| ^ 1. Since ~ e"'^'^', kz ~ sgn(y)e''l^l — 1, we find that the 
Schrodinger potential becomes 

V,{z)^—-—^ (A|y|»l). (3.18) 

(Mz\ + 1) 

^For a special case where A;A~^ is an integer, we can express the Schrodinger potential in terms of z explicitly. We show 
this in Appendix B 
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Figure 4: The Schrodinger potential for the scalar perturbation in terms of the coordinate system y. We choose 
j/o = 0, and k = 0.25. The solid line has A = 3 and the broken line has A = 1. 



This happens to be the same potential as that in the RandaU-Sundrum single wall model [H], in spite of 
different spacetime dimensions. The wave functions of the continuum massive modes for this potential 
are known to be expressed as linear combinations of Bessel functions at the region far from the wall [3] . 

3.3 The active scalar perturbation 

Next we study the perturbation of the active scalar field (^r. Notice that the fluctuation i^r around the 
active scalar field background 0r can be reduced to the trace part B of the metric perturbation through 
Ea. (|3.4(l . Therefore we mainly concentrate on the trace (scalar) part of the metric perturbation B in 
what follows. The linearized equation which contains only B can be derived by combining Ea. H3.3|l and 
()3.4|1 and using the background field equation : 



-2A 



□ 



D-l 



dl+[{D- 3)i -2^\dy + 2[D-i)[A- AX^ 



B^O. 



(3.19) 



In order to transform this into the Schrodinger form, we change the coordinate from y to z and redefine 

B. Substituting this into Ea. (|3.19|l . we find the Schodinger type equation 



the field as B 



-An ~1 



for the scalar perturbation; 

HeB = [-dl + Ve(z)] B = Ud-iB 
where the Schrodinger potential Ve{z) is defined by 

Ve(z)^-f^ + 2(^V + (^-4) 



D 



A 



(3.20) 



(3.21) 



Similarly to the tensor perturbation, the inner-product for the scalar perturbation B should be defined 
in terms of the conformal coordinate z and the redefined field B. 
Plugging our solution (|2.26() into this, we find 



Ve = [coshA(y - 2/0 ) 



A^ + M 1 



1 



cosh A(?; - yo) 



(3.22) 



We stress that Vg can be expressed in terms of y, but not in terms of z, since it is generally difficult to 
solve dz = e~^dy. This potential Ve has the following properties : i) it is positive definite, ii) it vanishes 
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asymptotically at infinity, and iii) the height of Ve is of order as shown in Fig^ From i) , it follows 
that there are no tachyonic modes since the wave function of such modes will necessarily diverge either at 
y = cx)ory = — cx). Therefore we can conclude that the background configuration H2.26|l is stable under 
the active scalar perturbation. From ii), it follows that the spectrum of the massive modes is continuous 
starting from zero. From iii), the potential diverges at any finite point y in the thin wall limit (A oo). 

Though we can not find the exact solutions for the massive KK modes, zero modes can be found by 
rewriting the Hamiltonian (|3.20() into SQM form as follows : 



log 



iAA' 



log 



0R, 



(3.23) 



To show this, we use the identity 0^ (A'" - 2A'A") = 2(/)^ [A" - A'^). Similarly to the tensor perturba- 
tion there are two zero modes of Tie ' 

B^°\x,z) = m\x) 4^e-^^ + B(-)(x) ^e-^^ [ dz ^e^D-2)A^ (3 34) 

Both zero modes are unphysical by the following reasons. The first term is unphysical, in the sense that 
this is eliminated by a gauge transformation preserving the Newton gauge (I3.1|l : 

p-{D-3)A „2A r 

6 = B^+Hx f-^-^, = -B(+)(,)^ j dy e-(^-i)^, (3.25) 

where is an infinitesimal coordinate transformation parameters. The transformation law is given 
in Appendix A. The second term is unphysical since it diverges at infinity and is not normalizable as 
illustrated in Fig[3 

Next we turn to analysis for the massive KK modes. As we have mentioned above, massive modes are 
continuous from zero. Similarly to the tensor perturbation, properties of mode functions can be examined 
by analyzing the behavior of the potential in the region far from the wall. In the region where lyjA 3> 1, 
the Schrodinger potential H3.22|l becomes : 

V,(z) c -^^±^. (3.26) 

This potential is very similar to the Schrodinger potential (|3.18|) for the tensor perturbation. Therefore 
all the massive modes are given by a linear combination of Bessel functions asymptotically at \z\ 00. 
Although these two Schrodinger potentials H3.18|l and H3.26|l have the same z dependence asymptotically 
|t/|A ^ 1, their behaviors in the thin wall limit are very different. The potential (|3.18(l depends only on 
k (fixed in the thin-wall limit), but not on A. On the other hand, the potential H3.26|l is proportional 
to polynomials in A. Therefore, the latter diverges in thin wall limit whereas the former is finite. This 
can be understood as follows. The perturbation of the trace part of the metric B is related to the active 
scalar field perturbation (^r through Ea. H3.4|l . Since all the massive KK modes associated with the active 
scalar field become infinitely heavy in the thin wall limit, the massive KK modes for the perturbation 
of the trace part of the metric freeze simultaneously. In this limit only the tensor perturbations remain 
which correspond to the known modes of the RS model® . 

^Generically speaking, the fluctuations of the inert scalar field <fii can be an exception depending on the potential, 
although the inert scalar {pi in our model is also frozen in the thin wall limit. 
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A' D-2 

' ■ : zero modes B'^'^\x)—-e 2—^ of the fluctuation of the trace part of the metric B in Ea. (|3.24|l 
can be translated into the perturbation of the active scalar field ipn by means of Eq. (|3.4|) : 

(x, y) = 0^+^ (x) ^Re-(^'3)A ^(+) (^^^Qy (3 27) 

where ip^\x) = [D — 2)b\i^\x) /2. In weak gravity limit [n 0), reduces to a constant. Then we 
find that this zero mode is localized on the wall and that it corresponds to the Nambu-Goldstone boson 
corresponding to the spontaneously broken translational invariance. 

3.4 Analysis for the perturbation about (f)\ 

In our tree level solution the imaginary part of the scalar field vanishes identically and does not 
contribute to the energy momentum tensor. Therefore it does not affect the spacetime geometry. We shall 
call scalar fields with no nontrivial field configuration as inert field. In the linear order of perturbations 
we found that the fiuctuation Lpi of this inert field decouples from any other fiuctuations, as shown in 
Ea.lTm. 

In order to find the spectrum of ipi, we first bring Ea. (j3.6(l into a Schrodinger form by changing the 
coordinate from y to z and redefining the field (pi = e^~^ipi. Then we obtain 



Hm = [-dl + Vi(z)] Do-m, Vi(z) = Vtiz) + e'^i^ 



(3.28) 



where Vt{z) is the potential for transverse traceless part of the metric defined in Ea. H3.8() . To obtain more 
concrete informations on the spectrum, we need to examine properties of each model. For our model we 
find 

-2M ^ ^' + ^(^ + 2cos^|0H)-^sm^f^H (3.29) 

We shall discuss generic property of this inert scalar for the non-BPS background in Sec^ 

If we choose the BPS solution as our background, we can rewrite the potential by using the BPS 
equations (|2.25(l 

1^ = + 2kA - 2A ~ 2A\ (3.30) 

Then the Schrodinger potential Vi takes the form : 

Vi = A'^ - A" + c^^ (A2 + 2fcA) . (3.31) 

We illustrate V/ in terms of y in FigEl For vanishing gravitational coupling k — > 0, Vi reduces to a 
constant A^, which agrees with the model of global SUSY in Ref. ^24. On the other hand, the potential 
Vi acquires regions of negative values when k becomes large. In the case of the BPS background, we 
can show that there are no tachyonic modes in this inert scalar sector with the aid of the SQM. Let us 
introduce a supercharge as Qi = —dz — A' and QI = dz — A'. Then the Hamiltonian Hi can be rewritten 



as 



Hi = qIQi + e^^ (A^ + 2fcA) . (3.32) 

The first term is a nonnegative definite Hermitian operator and the second term is never negative. 
Therefore, we can conclude that eigenvalues of Tii are always nonnegative and there are no tachyonic 
mode. 
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V 




Figure 5: The Schrodinger potential for the perturbation of the imaginary part of the scalar field in terms of the 
coordinate system y. We choose yo = and A = 1. The bold line has = 1, the solid line has k = 0.25 and the 
broken line has k — 0.05. 



4 Stability of Non-BPS multi- Walls 

For non-BPS solutions, the positivity of the energy of the fluctuation and the associated stabihty is entirely 
nontrivial. In the limit of vanishing gravitational coupling, however, our supergravity model reduces to a 
global SUSY model that has been shown to be stable [21] • Since the mass gap in the global SUSY model 
should not disappear even if we switch on the gravitational coupling infinitesimally, the massive scalar 
fluctuations in the global SUSY model should remain massive at least for small gravitational coupling. 
On the other hand, we need to watch out a possible new tachyonic instability associated with the metric 
fluctuations. 

As for the transverse traceless part of the metric, we have already shown that there are two possible 
zero mode candidates hl^'-'^\x) e^(^) and hl^'-~\x) e^(^) / dz e'^-^C^) in Ea. (|XT7| . In the non-BPS 
solution, we no longer need to worry about the normalizability of the wave function. Instead, we need 
to satisfy the boundary condition imposed by the presence of the boundary cosmological constants. To 
impose the boundary condition, we have to use coordinate system which is more more appropriate to 
specify the position of the boundary. This is achieved by going to the Gaussian normal coordinates [201 • 
We have been using the Newton gauge to study the mass spectrum in the Shrodinger type equation. We 
can follow the argument in Ref.pO] to obtain general coordinate transformations from the Newton 
gauge to the Gaussian normal gauge : 



1 



U^,y)^^ I dy'h22+£}^'{x), (4.1) 







1 



U^.v) = -^f dy'e-'^ r dy"h,,^,~^i^J r dy' e-''^ + £,f\x) , (4.2) 

^ .ly{±) Jy{±) Jy{±) 



where ^2^'',^^^'' depend on x only. Boundary conditions in the Newton gauge are found to be 



+2e— er;.. ' , =0' (4-3) 



(dy + A) h22 + 2ief' 1 = 0. (4.4) 
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We find that the former mode ft.JJ (x) e"*^^-* receives no constraint from the boundary condition, and 
is still a physical massless mode localized on the wall which should be regarded as the graviton in the 
effective theory. The latter mode hjjj^ •* (x) e'^'^) J dz e~^"*(^^ is constrained by the boundary condition 
()4.3|l . Since the constraint H4.3|l and (|4.4(l relate ^j^"* to h\ii} and 'h22{x), this mode should be classified as a 
scalar type perturbation. On the other hand, fi'^^i) can be gauged away through the gauge transformation 
(|3.25|) . which preserves the Newton gauge^. Therefore this mode becomes unphysical in the presence of 
the bulk scalar field like in our model. Since ^2 corresponds to the physical distance between the two 
branes |27j . our system automatically incorporates the stabilization mechanism without an additional 
bulk scalar fields. In the thin wall limit, the wall scalar field freezes out and the fluctuation of the scalar 
field ceases to be related to the scalar perturbation of the metric, so that h\j^i) can no longer be gauged 
away. As a result of restoration of h!-^J , any distance between two walls is admitted as classical solutions 
and the model becomes meta-stable. 

To evaluate the mass spectrum of massive modes, we use the small width approximation. Then the 
asymptotic behavior of the potential gives the wave functions expressed by means of the Bessel functions 
as in the Randall-Sundrum model. These eigenvalue spectrum is approximately equally spaced just like 
the plane wave solutions. These (almost) continuum modes should give the corrections to the effects of 
localized graviton, similarly to the Randall-Sundrum model. Therefore we obtain a massless graviton 
localized on the wall and a tower of massive KK modes for transverse traceless part of the metric. 

Since the active scalar field (^r exhibits a mass gap without any tachyon, we expect that there should 
be no tachyonic instability at least for small enough gravitational coupling. We have found in Ea. (|3.4ll 
that the active scalar field Lpii can be reduced to the trace part of the metric B. This implies that there 
should be no tachyon in the transverse traceless mode as well at least for small gravitational coupling, 
since both degrees of freedom represent one and the same dynamical degree of freedom. In fact, we have 
observed that the potential Ve defined in Ea. H3.21|l is everywhere positive and has no tachyon. 

There are two possible zero mode candidates as given in Ea. l|3.24|l . However, both of them are 
unphysical by the following reason. The first one can be eliminated by a gauge transformation. The 
second one has now no problem of normalization, since the extra dimension is now a finite interval, but 
it cannot satisfy the boundary conditions ^H]- 

To evaluate the mass of the lightest scalar particle, which is usually called radion, we use the thin-wall 
approximation where the wall width is assumed to be small compared to the radius of compactification R. 
To make this approximation, we separate the potential (|3.21|) for the trace part of the metric fluctuation 
into unperturbed and perturbation as 

Ve(.)^V(°)(.) + vW(^), (4.5) 



Vi°n^) = '/'ke'^^ , V«(z) = §-|2.^(</>^,)^ (4.6) 

where we considered in Z3-dimensions instead of 4 dimensions. The zcro-th order eigenfunction for the 
lowest eigenvalue is found to be 

B^'\z)^ , (4.7) 

e~2-^(/)J^ 



For simplicity, we have assumed Z2 parity of the metric perturbation to be even. 
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with the vanishing eigenvalue (mQ*^^)^ = and a normaUzation factor M . The first order eigenfunction is 
given by 

B^'\^) = f dz'e^^-^^^c^l f dz"e-(^^4)^<^;,-^ {yi'\z) - (m^)^) , (4.8) 

where the first correction to the mass squared eigenvalue is denoted as (toq^^)^. By using Ea. (|3.4() . the 
trace part of the metric can be transformed into the active scalar fluctuation (^r as 

where the position of the wall is at ?/ = and the boundary with the negative cosmological constant is 
y(+)- To satisfy the correct boundary conditions, we have to require that this first order eigenfunction 
(/5q^^ should vanish at the boundary . This determines the first order mass squared as 

ry(+) 



dy e~ 

«0^-2.^§^A___, (4.10) 

dy 



Taking 13 = 4 as in our model, and applying to the case of the first example of non-BPS background 
in Ea. p.H^II with the symmetry Z2 x Z2, we should identify = ■nR/2 and obtain 

« 8A2e-(i+"')^^^ (^1 + ^a^^ 2"^ ^ (4-11) 

where AT^ is the tension (energy density) of the wall. For the other background solution with the Z2 x Z2 
symmetry, we should identify the boundary with the negative cosmological constant as y(+) — 7ri?, and 
obtain the same result. It is interesting to note that the mass scale is given by the inverse wall width 
A, and that it becomes exponentially light as a function of the distance between the walls, even though 
the radion mass receives a complicated gravitational corrections. It is appropriate to fix the wall tension 
and the gravitational coupling in taking the small width limit Ai? — )■ 00. Then we obtain a simple mass 
formula in the limit ^ 

ml w SA^e-'^'^-^ ^ 0. (4.12) 

This characteristic feature of lightest massive scalar fluctuation is precisely the same as the global SUSY 
case Pni- The lightest massive mode in that case results from the fact that two walls have no commu- 
nication when they are far apart, and the translation zero modes of each wall becomes massless as the 
separation between walls goes to infinity. 

The mass spectrum of the inert scalar fluctuations Lpi is determined by the Schrordinger form of the 
eigenvalue problem (|3.28|l with the potential Vi{z). The potential has the same term as the transverse 
traceless mode Vt{z) with an additional term in Ea. H3.29(l which is nonnegative definite provided 

the gravitational coupling is not too strong k < ^ 

^ ^ (1-— 3-1 + 2^- >0. (4.13) 



9 J \ 9 



*This mass squared is factor two larger compared to the value of lightest massive scalar in the global SUSY model lTSl . 
We have not understood this discrepancy. 
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Therefore inert scalar does not produce any additional tachyonic instability. 



5 Fermions 



In the previous two sections, we focused on the stability of BPS and Non-BPS wall configurations and 
studied its fluctuations. In this section, we turn to the fcrmionic part of the model and study its 
fluctuation. We shall consider only the BPS solutions for simplicity, since it allows massless gravitino, 
whereas the non-BPS solutions do not. The part of the Lagrangian H2.1|l quadratic in fermion fields (with 
arbitrary powers of bosons) can be rewritten as 



— 1 ^quadratic 
^ fermion 



+e '"'"'' e 



V2 



-K^e'2 



V2 



+e-^^{D^Pxa^'4'^ + D^.P*x&tt^f.) + ^{D^PxaH2 + D^,P*x^H,2) 



V2' 



-{V^D^PX^ +V^.D^.P*X^) 



(5.1) 



The terms in the fourth line is quadratic in gravitino without any derivatives, which can be regarded 
as mass terms for gravitino. We find that they are Z2 odd under (j) — s- —(p. In this respect, our model 
provides an explicit realization of the condition to have a smooth limit of vanishing width of the wall 
|21| and in agreement with one version of the five-dimensional supergravity on the orbifold [TP . For our 
modified superpotential P given in Ea. H2.19|) . D^P and V^fjD^P are of the form: 

A3 „ a2 

r 



D^P 
V,,D^P 



A" g 
— cos— ( 
9 A 



2k^((^* - 4>)d^P~ I K- 

A2 



,.,2\ A^ . 9 



2k^ 



A2 g 
I — cos— <; 
9 A' 



(5.2) 
(5.3) 
(5.4) 



5.1 Gravitino 



In this subsection, we will explore a massless gravitino which is a superpartner of the massless localized 
graviton under the SUGRA transformation with the conserved Killing spinor (|2.27|) . Before studying 
equations of motion for gravitino, we will supertransform the wave function of the localized massless 
graviton to find conditions that the physical gravitino should satisfy. Let us focus on SUGRA transfor- 
mation law for vierbein in Ea. (|2.4(l . 



Seem- = in (Ccr-V'm + Cc^-V-m) • 

The preserved SUSY along the Killing spinor C{K) in Ea. (|2.27|l with 6* = 7r/2 is given by 



CiK) = -^Ca(if)a2""-^e■^/2[e2,-el], 



UK) 



.A/2f 



-ei,-e2\ 



(5.5) 

(5.6) 
(5.7) 
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Denoting the fluctuations hmn of the metric around the background spacetime metric gjjf^'^^'""""^ = 
diag{e^^rif^^,l) as gmn = + ^mn, the following Unearized 3D SUGRA transformations with 

the KiUing spinor C(^^) are obtained for the metric fluctuations 6hmn = S{em^na) = Sem-ena + emaScn- 

^C(.K)h^i, = iKe^CiK){cr^^|^^ - ia-a■^^p,, + a^-ip^ - ia-a^-i/jf^), (5.8) 
Sc(K)h2f, = iKC(iv:) |(T2V'At + + - »o-%^'(/;2)e'^| , (5.9) 

Sc(K)h22 = 2tKC{K){a2^p2+i^2)- (5.10) 

In sect 13. 21 we have imposed the gauge fixing condition (Newton Gauge) for graviton 

h22 - -^e-^^7f'h,,, = ^h, (5.11) 

/t2M = 0. (5.12) 

We can algebraically decompose /i^i/ into traceless part e^^hj^^ and trace part h. We have found that 
the localized graviton zero mode is contained in the traceless part 

77^'' V = 0- (5-13) 
Equation of motion shows that the localized graviton zero mode also satisfies the transverse condition: 

ri^^'dxh^. - 0. (5.14) 

The matter fermion of course do not have the graviton zero mode : = 0. 

It is useful to decompose Weyl spinors in four dimensions into two 2-component Majorana spinors in 
three dimensions. For instance gravitinos ipm are decomposed into two 2-component Majorana spinor- 
vectors ipm^ and ipln' (real and imaginary part of the Weyl spinor-vector) 

^ V'm. - ^'^L^t = -i^L^iP", (5.15) 

= + ^'^L^i = ^^LV^L'^"- (5.16) 

Similarly to the traceless and trace part decomposition of graviton (symmetric tensor), gravitino (vector- 
spinor) can also be algebraically decomposed into its traceless part i/ij and trace part "0 as 

V'm = V-J - ^^T^Vi, ^^i^l = o, S^V'm = ^- (5.17) 

Let us make a SUGRA transformations of the physical state conditions (|5.11(l - (|5.14|l for gravitons 
with the conserved KiUing spinor C(^''^) in Ea. l|5.7|l . The SUGRA transformations with ({K) of Eas. (|5.11|l . 
(jFm (|FT^ give 

= 6^^K) (^h22 + = 2tKC{K) (^as^'^ + ^e-^a/i^^f ^ , (5.18) 
= (5c(i^)/i2M = z<(a2^},'^+e%Vif ) . (5.19) 
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= <5c(K)77^''V = 2z«;e-^C'^ii^(,2). (5.20) 

These result suggest the most natural gauge fixing condition for local gauge SUGRA transformations 

= 0, (5.21) 

which can always be chosen. Then, the above gauge fixing conditions (|5.18I) - (|5.2U|I are translated as 
■0^^' — and the traceless condition for i]/^'^^ 

a^V^^ = 0. (5.22) 

Therefore we expect^ that the localized massless gravitino should be contained in the traceless part of 
-0^,^^ . The SUGRA transformation of the remaining condition (|5.14|l gives the transverse condition for 

(2) 

the iI^Il ■ Similarly to the graviton case, the localized gravitino should not have matter component 

X = 0. (5.23) 

Let us now examine the equations of motion for gravitino ip^^ coupled with the matter fermion 
which are obtained by varying the action H5.1(l . If we impose the conditions H5.21() , (|5.22|) , and H5.23|l on 
the gravitino equations of motion, we obtain 



(5.24) 



where we have used the BPS equation H2.15|l for background fields. Possible zero mode should give a 
vanishing eigenvalue for the operator in the parenthesis : 

= _eM2p^^^ J^Q^ + 1 + ^ _,^A..^2 (5 25) 

where e ^^'^^ Up = i{at^a-a—— rjt^a-) is used in the second equality. In terms of the 2-component Majorana 
spinors H5.15|l and (|5.16(l . we obtain 

(d2 + ^i) V^W^^ + (-d2 + ^i) ^(2)^^ = 0. (5.26) 



Since ipj^^ = 0, we obtain 



-d2 + ii) V?)^^ = 0. (5.27) 



Therefore, we find the gravitino zero mode in the transverse traceless part of the 2-component Majorana 
vector-spinor ■0^ with the wave function 

f;!'^^{y) = e^. (5.28) 

^One should have in mind that it is desirable to choose a gauge fixing condition for SUGRA transformations to become 
a supertransformation under C,(K) of the gauge fixing condition for general coordinate transformations. However, it may 
not be logically mandatory. 
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Now we see that the locahzed massless gravitino wave function is in precise agreement with that 
expected from the preserved SUGRA transformation with the Killing spinor : 



(5.29) 



Since the wave function of the graviton and the Killing spinor are e^— ~ e"* and ({K) = e2 , we find 



5.2 Matter Fermion 

In this subsection, we study the fluctuation of matter fermion. By varying the Lagrangian (|5.1|) with 
respect to x, we obtain the equation of motion for matter fermion x- Using the gauge choice ip2 = '4'2 — 
to the equation of motion, we find 



= -ie^a^d^x-^iA^2X-^^-d2X-e'^'^''^'DrDrP*X 



(5.30) 



The second line gives the mixing term between the trace part of gravitino ip = <7'^'ipf_i and the matter 
fermion x- Using the BPS equation H2.18|l . the mixing term can be rewritten as 



^--v-M ^ 



^a!ia^^ - ^^a^i^f, = -^<t>^- (^m " = ~— -^^-^1'^ (5-31) 



where we used the 2-component Majorana spinor notation defined in Eas. l|5.15(l . (|5.16|l . Since the mixing 
occurs only with it is also useful to decompose matter fermions into 2-component Majorana spinors, 
similarly to Eas. H5.15|) . H5.16|l . Then the matter equation of motion is decomposed into two parts with 
opposite transformation property under the charge conjugation cr^ 

3 



d2 + -A~e'^^''^V^,D^,P* 



X 



(1) 



(5.32) 



-le a- 



d2 + -A + e'^'^''^V^,D^,P* 



It is now clear that we have a zero mode consisting of purely x'^'' : 
The zero mode wave function for matter fermion is given by 



d2 + -A-e-'^ ^V^^D^.P* 



= 0, 



whose solution is given by 



x[,^)~e-3^/2exp 



dy e" 



X 



'V^D^P 



(5.33) 



(5.34) 



(5.35) 



(5.36) 



Using (|5.4() . the integral in H5.36|l in BPS case reads 

dy e'^^'^V^D^P = -{k^ + j^)^ J dy sin|0 = -(1 + «:'^)A / dy tanh(Ay) 

= -(1 + ^2 — )log(cosh(A2/)). 
9 



(5.37) 
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In the second equaUty, P — e '^^'^^ ^sin{g(j)/ A) in Ea. H2.19|) is substituted and cj) — (j)* is taken into 



g 

account. In the last equality, the BPS solution 4> = j {2taii-'^e^'^y-y°^ - f ) in Ea. lfT^ with n = is 

considered. Then the zero mode wave function of matter fermion x^^^ is given by 

X^,'^ ~ e-3^/2[^Qgl^(^y)p[iW^] _ ^5 33^ 

In the weak gravity limit, the zero mode of matter fermion reduces to the Nambu-Goldstone fermion 
associated with the spontaneously broken SUSY ^3] 

Xo^^ V «^0. (5.39) 

cosh(Ay) 

As expected in the global SUSY limit, the wave function is localized at the wall where two out of four 
SUSY are broken. Let us note, however, that this zero mode of matter fermion should be unphysical 
except at K = limit. For any finite values of k, it should be possible to gauge away this zero mode, 
precisely analogously to the zero modes in Eg. (|3. 24(1 in the matter scalar sector in sect 13. 31 In fact 
we can see that the A dependence (warp factor) of the zero modes of active scalar 5^+^ (x) and the matter 
fermion x'oH^) agrees with the surviving SUSY transformation generated by the Killing spinor C(^)) 
and will form a supermultiplet under the surviving SUGRA, since (f>^\y) ~ s~^iXo^\y) e~^^^'^ and 

C(y) - 

<5c0(x,j/)r = V2ax,y)/'\x,y). (5.40) 

On the other hand, x^"^"^ should contain another Nambu-Goldstone fermion corresponding to the SUSY 
charges broken by the negative tension brane, if we consider non-BPS multi-wall configurations. Noting 
that the mixing term is suppressed by the Planck scale Mp, the zero mode equation of motion for x'-^'' 
in weak gravity limit k — s- is given by 

= d2X^^'> +e'^'^''^V^D^Px^^\ (5.41) 
^ 82/^^ - Asin|-0 x^^^ = d2X.^^^ - Atanh(Ay)x'^\ (k ^ 0) (5.42) 

where BPS solution = (2tan^^e'^(^^^''^ — ^) is substituted in the last equality, thus the zero mode 
wave function becomes 

xl,'^ ^cosh(A2/) (k^O), (5.43) 

which is not normalizable and hence unphysical even in the limit of k — > 0. We know from the exact 
solution of the non-BPS two- wall solution|15|. that this wave function results when taking the limit of 
large radius to obtain the BPS solution. In that limit, SUSY broken on the second wall at y = irR 
is restored and the corresponding Nambu-Goldstone fermion, which was localized on the second brane, 
becomes non- normalizable and unphysical. This is precisely our zero mode wave function Xq ■ 
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A Appendix A 



In this appendix we show the gauge fixing to the Newton gauge. The most general fluctuation around 
the background metric (|2.6I) takes the form : 

ds^ = e^'^{Tj^, + hl^ + 2i^f,,B)dx''dx'' + 2ff,dx''dy+il-2C)dy\ (A.l) 

where the trace part of the (/i, v) component of fluctuations is denoted as 27y^jyi3 and the traceless part 
is denoted as /ijj^. denotes the fluctuation of (/i, 2) component and — 2C denotes the fluctuation of 
(2, 2) component. 

After a tedious calculation we find the linearized Ricci tensors : 



<J - 1^'^ (e-^-^DD-i + + {D~ l)Ady + 2{D - 1)^^ + 2A) - hl^^,^ 
+V^.^e^^ (e'^^DD-i +dl + 2{D - 1) (Ady + i^) +2A^B + {D- 3)B 
-V^.Afp" - [dy + {D~ 3)i) + rif,,e^^ [Ady + 2A + 2{D - l)A^'^ C - C^,,, (A.2) 



<i = -\ (-c-^-^Dz,-! - 2{D l)A^ 2A) U ^c-'^f,.,^" 

HD - 2)AC^ + {D^ 2)dyB^^ - \dyhl/, (A.3) 

^2' - - 1) [dl + 2Ady) B - c-^^dyfp^P - -{D- C, (A.4) 

where we define B^^ = d^B, hj^^ ,-^ = ^ (^-^^.j. + ^i^p,^)' ^p'' ^ ''f^^>^Ip and Ud-i = rf^dpdx. We also 
find the linearized energy momentum tensor as follows : 



f (1) ^ ^ 

D-2 



(A.5) 

f^^^ = 2<^R^R^^ + -^^VuU, (A.6) 

t4^) = 4<^h5,^r, + ^^^R - ^^rC, (A.7) 

where ip^i is the fluctuation around the background active scalar field i^r. Notice that the fluctuation 
ipi about the background configuration for the imaginary part (pi decouples from any other fields in 
linear order of the fluctuations. We can obtain the linearized Einstein equations by plugging these into 

_„.2f(l) 

The above results are the most general in the sense that we do not fix any gauge for the fluctuations. 
As a next step, we wish to fix the gauge that simplifies the linearized equations. The "Newton" gauge is 
known as a candidate of such a gauge |3()l 02] . The gauge transformation laws for the fluctuations are of 
the form : 

Shi, - -^(p..) + -^^Vpuip", SB = -i6 - ;^^Cp'^ 

(5^ = -e2^a.„4-6.p, SC = dy^2, <^^r--0r6, (A.8) 

where is an infinitesimal coordinate transformation parameter Sxm = £,m and = e^^'^^^. Using 
these four gauge freedom, we fix = and {D — 3)B = C. The residual gauge transformation should 
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satisfy 

dyi^ + e-2^C2,M - 0, (9, + {D- 3)i) 6 = -^^ip '- (A.9) 
In this gauge the hnearized Einstein equations take the form : 

ie^-^ {e-'^Dn-^ + + (D - l)Ady + 2{D 1)A' + 2A) - h^^^^y" + ^rj.^hj,'^'' 

^ n^^^VA. (A.IO) 



£1-2 



e"^ { e-^^Uo-i + dl + {iD - b)Ady + 2{D - 2) ({D - l)A^ + ^)) ^ - -^—fix'"''' 



^_-.V-^2^H,5 + -^^Hj, (A.11) 
(p - 2) (a^ + (I? - 3)i) B,^ - is./iT^-" = -2k2^r(^r,^, (A.12) 

where the Ea. l|A.10|l is the tracclcss part of the (/i, v) component whereas the Ea. ljA.lip is the trace part 
of it. Denoting D{x, y) ^ {D - 2) (dy + {D - 3)A) B + 2k20r(^r, the Ea. (|A.12|) is rewritten as 



Summing the background Einstein equation Ea. (|A.ll() muhiphed by {D — 3)e~^'^ and Ea. l|A.13p gives 
(dy + {D- D{x, y) = 2(I)~ V "'^^^'^''''^' ^"""^ 

= e-(^-i)^(^) + ^^^e-(^-i)^| dy e^^-^)^/.^^^-^ (A.15) 



where -E(a:;) is an arbitrary function of x. At this stage, the hnearized Einstein equations are Ea. ljA.10|l . 
(Enj, and (f05)l . 

Next we attempt to ehminate the longitudinal mode of /ij^ by using the residual gauge freedom. That 
is, we wish to set = hj^p''' — 0. For that purpose we first derive the equations of motion for u^. Taking 
a divergence of Eq. (|A.10|I , we find 

(dy + {D- l)A) dyv, = ^e-2^^;,-V. (A.16) 

This equation can be solved as follows : i) taking divergence, we can determine Vp'^, ii) regarding the 
solution Vp''' as a source, we can determine w^. The gauge transformation law of takes the form 

Sv^ = -Do-ii^ - ^^ip,, ". (A. 17) 



We want to set = Vf^+Sv^j, by using the gauge transformation ljA.17p whose satisfies the condition ljA.9|l 
for the residual gauge transformation. Notice that the equations for Vf^ and dVf^ are identical second order 
differential equations since the gauge transformation law consistent with the gauge condition l|A.9|l does 
not change the form of the equation (|A.16|I . We can also verify this from the condition (|A.9|I straight- 
forwardly as follows. From Eg. lfO)! we find the identity e-'^'^DD-i^2 = (dy + {D ~ 3)i) 6- 
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Combining this and Kn. fOTIl . we find Svp^P = ^^^^^^p^e^-^ (^dy + (D - l)i) dy (^dy + {D - 3)i) $2 
and dySv^ = '^''0-3 ^ i^y + {D — 3)Aj Hence, (5w^t satisfies just the same equation as Ea. (|A.16|l : 



[dy + {D- l)i) dySv, = ^-^e-'^Sv,^^,,. (A.18) 

Therefore, can be ehminatcd in the gauge, if we can set at a given y = yo surface 

Vp, = -Sv^, dyVf^, = -dydv^,. (A. 19) 

To clear matters, we introduce new functions Ap(x) = ^^(x, yo), ^^{x) = (a;, yo), r(x) = 1^2(2^, yo), 

A{x) = {dyS,2) {x, yo) which are defined at y = yo surface. In terms of these functions Ea. l|A.19p can be 
rewritten as 



□i3-iA,, + ^-|Ap,,,^'' = ^^, D^^iS^ + = 6^, (A.20) 



where A^{x) = v^{x,yo) and Bf^{x) — {dyV ^) {x , yo) ■ Similarly, the gauge condition ljA.9|l at y = yo 
surface can be rewritten as 

+ = 0, A + (z?-3)ir = -^5|Ap'^. (A.21) 

Ap and can be determined similarly to the Ea. l)A.16|l . Next, we determine T from the first equation of 
Ea. (|A.21|l . However, this equation does not necessarily have a solution for a general function S^. To see 
this in detail, plug this into the second equation of Ea. l|A.20|l and we obtain UId^iT = — ^'^Ak 



^ 2L>-5 ^' 

This equation can be solved if and only if is expressed as a gradient of some function. In our 
case we obtain = "id^D from Ea. ljA.14p . Hence, a solution F of the first equation of Eg. (|A. 21(1 
exists. At the end, A is determined from the second equation of Ea. (|A.21|) . In this gauge, we obtain 
D{x, y) — i?e^'^~^^"^ where i? is a constant from Ea. (|A.14|) and (|A.15|) . We set E — since we require 
that the fluctuations B and </3r should vanish at infinity |a;| — > cxd on the wall. Thus we established our 
gauge choice (Newton gauge) and the constraints for the residual gauge transformations are l(A.9(l and 

□d-iCm + §3|ep'% = 0. (A.22) 

B Appendix B 

For a special case where fcA^^ is an integer"'^" , we can express the Schrodinger potential in terms of z 
explicitly. As an illustrative example, let us take fc = A, where we find (putting yo ~ 0) e"^ — (cosh /cy)"^, 
z — k^^ sinhfcy, A{z) — — - log (l + k^z^). The Schrodinger potential for the tensor perturbation takes 
the form : 

fc^(l-2fcV) 

There remains only one parameter controlling both the width of the wall and the magnitude of the 
gravitational coupling, similarly to r!,ef.|Hl|. Zero mode wave functions can also be expressed in terms of 
the z coordinate explicitly and are shown in Fig|Sl : 

/,TT(0)(^ ) ^ C^(+)(x) , ^ + hlJ^-\x) (B.2) 
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Figure 6: The Schrodinger potential for the tensor perturbation in terms of the coordinate system z. We choose 
fc = A = 1. The bold line denotes the potential, the solid line denotes the graviton fe^J'^' and the broken line denotes 

'TT(-) 
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Figure 7: The Schrodinger potential for the scalar perturbation in terms of the coordinate system z. We choose 
fc = A = 1. The bold line denotes the potential, the solid line denotes i3'+^ and the broken line denotes B^-^\ 



The Schrodinger potential Vg for the active scalar perturbation can be also expressed in terms of z : 

e (3 + 2fc2^2-) 



(B.3) 



(1 + ^22)^ 

and zero modes are of the form : 

i3(°)(a;,z) = zVTTfc2^ + i3^"'(a;) [l + zi-AiC^ kz) . (B.4) 

These are shown in FigEI 
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